A Renormalization Group For Treating 2D Coupled Arrays of Continuum ID Systems 



O 

o 

(N 



00 

h 



I 

G 

O 

o 



> 
o 



o 
o 



X 



Robert M. Konik 1 and Yury Adamov 2 

1 CMPMS Department, Brookhaven National Laboratory, Upton, NY 11973 
2 Department of Physics, Texas A&M University, College Station, TX 77843 

(Dated: June 23, 2007) 

We study the spectrum of two dimensional coupled arrays of continuum one-dimensional systems 
by wedding a density matrix renormalization group procedure to a renormalization group improved 
truncated spectrum approach. To illustrate the approach we study the spectrum of large arrays of 
coupled quantum Ising chains. We demonstrate explicitly that the method can treat the various 
regimes of chains, in particular the three dimensional Ising ordering transition the chains undergo 
as a function of interchain coupling. 
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The density matrix renormalization group (DMRG) [l[ 
is one of the primary theoretical tools for the quantitative 
description of low dimensional lattice models. For a wide 
range of one dimensional (ID) lattice models, DMRG 
can characterize the model's spectrum and correlation 
functions Q. While there have been notable recent ad- 
vancements 0,0], its use on 2D lattice models is more 
circumscribed |5|]. 

There exist several strategies to apply DMRG to 2D 
models. In the first, the 2D lattice is reduced to a ID lat- 
tice with long range interactions @. A second approach 
sees short chains treated as individual lattice sites, allow- 
ing the ID DMRG algorithm to be applied to a model 
with short ranged interactions directly Q. In a more 
sophisticated variant of this methodology, the DMRG is 
applied in a two stage process Q . The 2D system is first 
divided into a set of coupled ID chains and the DMRG 
is used to determine a low energy reduction thereof. For 
the second stage, the reduced chains, coupled together 
and treated as individual lattice sites in a ID lattice, are 
analyzed again using DMRG. In a final approach, the ID 
matrix product states underlying the DMRG algorithm 
Q are replaced by a higher dimensional generalization, 
projected entangled pair states 0. 

In this letter we present a distinct approach to apply- 
ing DMRG to 2D models. This approach trades upon a 
description of a 2D system as a mixture of continuum and 
discrete degrees of freedom. In particular, we approach 
2D systems as coupled arrays of continuum ID chains 
with truncated Hilbert spaces. This methodology offers 
several distinct advantages. It allows us to treat any 2D 
strongly correlated model provided it can be conceived as 
composed of continuum ID subunits. Furthermore, the 
approach affords superior finite size scaling. As a func- 
tion of the length, R, of the composite ID systems, finite 
size corrections behave exponentially. This implies that 
we can access, at the very least, the infinite volume limit 
in the dimension parallel to the chains. Finally, the trun- 
cation of the underlying ID Hilbert space dramatically 
lessens the numerical burden of the DMRG algorithm, 
while providing a natural means to perform a Wilsonian 



renormalization group (RG) improvement of any result- 
ing answer [Io| . 

The specific type of system that we propose to study 
takes the form, 



(1) 
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The ID continuum subunits of the array are governed 
by H} D which we insist must be either gapless (and so 
governed by a conformal field theory) or gapped but in- 
tcgrable Thus this method can study, for exam- 

ple, arrays of Luttinger liquids and a wide range of cou- 
pled ID Mott insulators. The subunits are coupled to- 
gether via the nearest neighbour coupling JOiOj, where 
Oi is an operator defined along the i-th continuum chain. 
This coupling should be relevant but can be of arbitrary 
strength. 

The analysis of the arrays proceeds in two conceptual 
steps. In the first step we follow Zamolodchikov's pio- 
neering numerical analysis of perturbed gapless ID con- 
tinuum theories an approach termed the truncated 
spectrum approach (TSA). We thus place the ID chains 
on a ring of circumference, R. Unlike DMRG applied to 
pure lattice models, periodic boundary conditions along 
the chains can be employed without issue. By working at 
finite R, we discretize the ID spectrum. This permits the 
states in the chains' spectrum to be ordered in energy, i.e. 
|1), |2), . . ., and then truncated at some finite cutoff, E c , 
leaving us with a finite number of states in the theory. 
With these alterations we nonetheless remain in an ex- 
cellent position to obtain information regarding the full 
theory in infinite volume. In Ref. [l2l ], a critical Ising 
chain in a magnetic field was studied. Choosing E c so 
that a mere 39 states were kept, infinite volume results 
were reproduced within an error of 2% (via diagonaliza- 
tion of a 39 x 39 matrix). It was this finding of excellent 
results at little numerical cost that motivated us to apply 
the TSA to more complicated situations [Io| . including, 
as here, arrays of ID chains. 

A part of the TSA's success is predicated on embed- 
ding non-perturbative information into the initial com- 
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TABLE I: Finite system DMRG algorithm adapted to the 
presence of a truncated spectrum. 

1. Form initial Hamiltonian, ff m _i, (and any other needed 
operators, O m -i) of system block, B m _i, of m-1 chains. 

2. Form Hamiltonian of superblock, fl m _i • •B m -i, of 2m 
chains, only keeping states whose energy, governed by 
H m -\ ® H 1D ® H 1D Cg) Hm-i, i.e. coupling between 
the blocks is absent, is less than E c . 

3. Find low-energy target state(s) of superblock. Form 
reduced density matrix of m-chain block, p m . 

4. Diagonalize p m , keeping sufficient eigenstates to 
obtain a truncation error of less than e c . 

5. Recast H m and O m in basis of kept eigenstates of p m , 
obtaining H m and O m . 

6. Diagonalize Hm, obtaining H m . Rewrite O m in terms of 
eigenstates of H* n . 

7. Form new superblock, B m • •B m , of length 2m+2. 

As determined by H m ® H 1D ® H 1D ® H m , only keep 
states whose energy is less than E c . 

8. Repeat 3-7 with m — » m + 1 until desired length of 
system is obtained. 

9. Perform finite volume sweeps until convergence. 



putation. In the context at hand this means that given 
any two states, |fc), |Z), we can exactly compute matrix 
elements of the form (k\Oi\l) by virtue of H\ D either be- 
ing conformal or integrable. We stress that computation 
of such matrix elements is always a practical possibil- 
ity either by exploiting the algebraic structures inherent 
in a conformal field theory or the form factor bootstrap 
approach in the case of a integrable field theory (l3j . 

Having prepared the chains by truncating their finite 
volume spectrum, we proceed to the actual analysis of 
coupled arrays. We perform the analysis using the finite 
volume DMRG algorithm. Here each chain with trun- 
cated spectrum is treated as if a single site. We however 
adapt the DMRG algorithm to take into account that we 
are working with a truncated spectrum. This algorithm 
is outlined in Table 1. 

Most steps of our finite system DMRG alg orithm are 
unchanged from the standard algorithm [l4 |. One pri- 
mary difference lies in the formation of the superblock 
Hilbert space and Hamiltonian (steps 2 and 7 in Table 
1). Instead of keeping all states in the Hilbert space, we 
insist that the state have an energy less than the trun- 
cation energy, E c , as governed by the Hamiltonian of 
the uncoupled blocks. Concomitantly in order to mean- 
ingfully associate an energy with an arbitrary state of a 
superblock, we must be able to assign an energy to the 
states of each block. 

This mandates that at each iteration of the DMRG 
we must not only recast the block Hamiltonian, H m , in 
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FIG. 1: a) A sketch of an array of Ising chains, each of length 
R; b) a zero temperature phase diagram of the coupled chains 
showing both the ordered and disordered regions. 



terms of the kept eigenstates of the reduced density ma- 
trix, p m (step 5), obtaining H m , but we must diagonalize 
H m (step 6). Working then with a basis of eigenstates 
of H m allows us to associate a definite energy to the su- 
perblock states formed in the next step (step 7). 

A notable operational feature of this DMRG is the rel- 
atively small number of states we need to keep from the 
diagonalization of the reduced density matrix in order to 
achieve a given truncation error. In the example of arrays 
of coupled quantum Ising chains that we consider below, 
the number of kept states needed to achieve truncation 
errors on the order of 5 x 10 -5 ranges from 10 — 40 deep 
in the ordered phase of the chains to 50 — 90 near the 
critical value of the interchain coupling where the chains 
order. 

The relatively small number of states needed to achieve 
a given accuracy mimics Zamolodchikov's original finding 
that a relatively small number of states was needed to 
describe a critical Ising chain in a magnetic field. But it 
also reflects the presence of a truncated spectrum in our 
approach. As was shown in Ref. fl5j |. the entanglement 
entropy that arises from dividing a 2D system into two 
behaves as the cutoff, E c . As the entanglement entropy 
is one measure of whether a DMRG-likc algorithm will 
be successful [H], ||| , we believe that our introduction of a 
cutoff into the problem is a key feature of our approach. 

A second key but related feature that marks our use of 
the DMRG algorithm is the use of RG improvement. For 
sufficiently large truncation energies, E c , the quantities of 
interest (whether it be the spectrum or some observable) 
satisfy a one- loop RG flow [10( : 



dAQ 
d\nE r 



= -gAQ. 



(2) 



Here AQ = Q{E C ) — Q{E C = oo) is the deviation of 
some quantity as a function of the truncation energy from 
its value at E c = +oo. The function g is a constant 
determinable from high energy perturbation theory and 
is related to the anomalous dimension of Q. When Q is 
an energy of some state, g = 1. Knowledge of this flow 
allows us to run the DMRG algorithm at several different 
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J - Interchain Coupling 



FIG. 2: The ground state and first excited state energy of a 
100 chain array as a function of interchain coupling obtained 
with both a random phase approximation (blue solid) and our 
DMRG methodology (red dashed). Here A = 1. 
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FIG. 3: The behavior of the ground state energy as a function 
of chain length, R, for 100 chains. J± = 0.275. In blue we see 
the analytic prediction while in red we see the results of the 
DMRG. 



truncation energies and then extrapolate the resulting 
flow to E c — +oo, so removing any residual effect of our 
use of a finite value for E c . 

Coupled Ising Chains: To demonstrate this method- 
ology we consider arrays of Ising chains (Fig. la) coupled 
via a nearest spin-spin perturbation: 

H = J2Hl DIsin9 ~JxJ2^h ( 3 ) 

i <ij> 

where the summations are over the chains in the array. 
The lattice form of the ID Ising model is 

h xd Ising = _j^tat +1 + (1 + g)af). (4) 

i 

Its continuum version is a Majorana fermion with gap, 
A = gj . We place the chain on a ring of circumference, 
R. The corresponding Hilbcrt space of the chain divides 
itself into two sectors, termed Neveu-Schwarz (NS) and 
Ramond (R). At T=0, the chain can be either ordered 
or disordered. In its disordered phase, A < 0, the NS/R 
sectors of the chain permit only even/odd numbers of free 
fermionic excitations. In the ordered phase, A > 0, the 
two sectors both permit only states with even numbers 
of fermions. The matrix elements of the spin operator, 
(i\a z \j), needed to carry out the DMRG procedure can be 
found in Ref. [loT ]. These matrix elements are only non- 
vanishing if the states \i) and \j) lie in different sectors of 
the theory. Because we couple the chains together with 
nearest neighbour spin bilinears, i.e. of of . 1; the Hilbcrt 
space of the full theory possesses two sectors. Any given 
state of the full theory has a tensor form C3i|fcj) where 
\ki) is a state on the i-th chain. These two sectors are 
distinguished by whether an even or odd number of the 
\ki) lie in the NS sector (or equivalently, the R sector) 
of the individual chains. That the Hilbert space possess 
different sectors is a generic feature of continuum models 
when placed in finite volume and is not particular to the 
Ising chains at hand. 



We now submit the DMRG algorithm to a number of 
tests. The first test that we apply relates to the be- 
haviour of the spectrum of the chains deep in their or- 
dered phase (A > - see the phase diagram in Figure 
lb). In this region of the phase diagram, we expect a 
chain RPA analysis to be accurate [l8j and so provide a 
baseline of comparison for our DMRG results. The chain 
RPA analysis amounts to treating the model 

Hurray = j £ ^ ^ + ( j + g)gf + ^ ^ (g) 

i 

where h RPA is chosen in the standard self-consistent fash- 
ion [3]. We analyze this model numerically using the 
truncated spectrum approach along the same lines as Ref. 
[l6| . However akin to the discussion surrounding Eqn. 
(2) and in Ref. [Io| . we perform an RG improvement of 
our numerical results. 

In Fig. 2 we plot the ground and first excited state en- 
ergy of a 100-chain array as a function of the interchain 
coupling as computed by both the DMRG and the chain 
RPA analysis. In order to optimize numerical perfor- 
mance, we employ different chains lengths for different 
values of J±. To obtain R = oo values for the gaps, 
A exc , of single excitations, we can use any finite value of 
R provided that we satisfy RA exc ^> 1. Because finite 
R corrections behave exponentially, 6rA bxc ~ e --RA ea;c ^ 
this constraint need only be satisfied loosely. In order 
to obtain a truncation error, e c of 1 x 10 -4 we needed to 
keep at most 18 states. Decreasing the truncation er- 
ror to 2.5 x 10~ 6 we needed to keep at most 34 states. 
Moreover by comparing our results with an analysis of 50 
coupled chains, we know that any finite size error related 
to studying only a finite number of chains is extremely 
small. The only significant uncertainty comes from the 
RG improvement - applying Eqn. (2) - and is the source 
of the error bars in Fig. 2. 

We see from Fig. 2 that the ground state energy as 
determined by the DMRG agrees exceptionally well with 
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FIG. 4: The behavior of the gap in the disordered phase as 
a function of J± for 60 coupled chains at R = 10 A. In red 
we plot the gap as determined for the truncation energy E c — 
7.8A, while in blue we plot the RG-improved curve. For both, 
the values of J± c and v are indicated. 

the chain RPA analysis even for values of the interchain 
coupling on the same order as the gap. The DMRG val- 
ues of the gap for the first excited state however see in- 
creasing deviations from the RPA values as the inter- 
chain coupling is increased. These results indicate that 
the truncated DMRG algorithm is operating as expected. 

A more significant test of the algorithm is whether it 
predicts correctly the finite R corrections. These correc- 
tions can be computed analytically, at least at leading 
order. According to Ref. [Ill, they take the form 

^ = -^EE/*^" fi£b(Ml) ' ( 6 ) 

b k y 

and can be thought of as arising from spontaneous emis- 
sion of excitations from the vacuum. Here eb{k y ,k x ) is 
the dispersion of some band b of excitations in the cou- 
pled chains. 8nE gs sees contributions from excitations 
with different (discrete) values of momenta, k y , trans- 
verse to the chains as well as excitations with different 
(continuum) values of the momenta, k x , parallel to the 
chains. While we cannot compute e exactly, we can esti- 
mate it from the values of the excited states as measured 
by the DMRG together with an RPA analysis [3]: 

e b {k y ,k x ) = (v 2 Fx k 2 x + A 2 exc + 2J ± Z (l-cos(k y ))) 1 / 2 . 

(7) 

Here Af, exc marks the lowest lying excitation in the 
band, while Z Q is the square of a matrix element of the 
spin operator on a single chain, Z a = 2 (0 of (0)|e&(A; x = 
0,k y = 0)) 2 . The latter quantity is estimated from the 
same analysis used in computing the RPA energies of Fig. 
1. 

In Fig. 3 we plot E gs {R) for a particular value of inter- 
chain coupling (Jj_ = 0.275) using the correction in Eqn. 
6 by taking into account the two lowest energy bands in 
the theory (blue dashed line) vs. the direct DMRG com- 
putation (red line). We see that for all but the smallest 



values of R — where our leading order analytic approxi- 
mation breaks down - the DMRG and expected analytic 
values agree well. For a point of comparison, we also 
plot the extrapolation of the energy from large values 
of R, where E gs scales linearly with system volume, i.e. 
E gs {R) = a(J±)RN 'chain (straight dashed orange curve). 
For values of R ~ Ar* the finite size corrections are 
exponentially suppressed, consistent with A& e xc = 5.1A 
for Jj_ = 0.275 (from Fig. 1). 

The final test we put to the truncated DMRG algo- 
rithm is the chains' ordering transition. Beginning with 
chains in their disordered state (A < 0) and coupling 
them together with increasing Jj_, they eventually order. 
This transition is in the same universality class as the 
3D Ising model. This implies that the gap, A exc , in the 
disordered phase should vanish as A exc ~ (Jj_ c — J±) v 
with v = .630 [3]. From our DMRG analysis (sec Fig. 
4), we find after RG improvement, J± c = .184 ± .0025, 
together with good agreement for the critical exponent, 
v = .622 ± .019. We see that RG improvement notably 
improves the results. Computing instead Jj_ c and v from 
the results obtained at the largest of the truncation en- 
ergies employed (E c = 7.8A), we obtain J c = .1880 and 
v = 0.650. 

This DMRG approach to 2D arrays of continuum sys- 
tems has a number of potentially valuable variations. We 
first note that we have managed to treat 2D arrays us- 
ing essentially a ID DMRG algorithm. Ref. |3J] has 
demonstrated that the DMRG algorithm, if used in a two 
stage process, can treat systems of one higher dimension. 
This implies that it may well be possible to study 3D 
arrays using our approach. We also note that the algo- 
rithm we have outlined in this letter may see substantial 
improvement if wedded to a numerical rcnormalization 
group (NRG) a la Wilson 0. At each DMRG itera- 
tion, one could perform a NRG akin to that described in 
Ref. (Io| to dramatically increase the truncation energy 
being employed and so (hopefully) dramatically improve 
the results of the procedure. 

To summarize, we have demonstrated that the ID 
DMRG algorithm can be directly applied to 2D arrays 
of ID continuum systems. In particular we have shown 
this algorithm can describe the behavior of large arrays of 
quantum Ising chains both in their ordered phase and in 
the vicinity of their order-disorder transition. We expect 
that this procedure will produce quantitatively accurate 
results on a wide variety of 2D systems in their infinite 
volume limit. 

RMK and YA acknowledge support from the US DOE 
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